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Abstract. 

The crystal structures and the physical (magnetic, electrical transport and 
thermodynamic) properties of the ternary compounds CeRhSi2 and Ce2Rh3Si5 
(orthorhombic CeNiSi2- and U2C03 Sis-type structures, respectively) were studied in 
wide ranges of temperature and magnetic field strength. The results revealed that 
both materials are valence fluctuating systems, in line with previous literature reports. 
Direct evidence for valence fluctuations was obtained by means of Ce Lm-edge x-ray 
absorption spectroscopy and Ce 3d core-level x-ray photoelectron spectroscopy. The 
experimental data were confronted with the results of ab-initio calculations of the 
electronic band structures in both compounds. 
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1. Introduction 

The system Ce-Rh-Si is known to be exceptionally rich in ternary phases, which scan a 
full spectrum of magnetic behavior related to the hybridization of cerium 4/ electronic 
states with d and p states of neighboring atoms. The most intensively studied has been 
so far the tetragonal compound CeRh2Si2 that is considered as an archetypal example 
of pressure-induced superconductivity emerging from magnetically ordered state \\\ . At 
ambient pressure, the 4/ electrons in CeRh2Si2 are well localized and give rise to an 
antiferromagnetic ordering setting in at Tni = 36 K that is followed by a change of 
the magnetic structure at Tn2 = 25 K [2]. Under high pressure one observes a gradual 
suppression of the Neel temperature down to absolute zero in the vicinity of the critical 
pressure Pc ~ 1 GPa |3]. Most interestingly, in the quantum critical region near the 
compound becomes superconducting below Tc = 0.5 K [3lll]. Simultaneously, the Fermi 
surface topology varies abruptly and the cyclotron mass measured in de Haas - van 
Alphen experiment increases rapidly, hence manifesting a change in character of the 4/ 
electrons, which become itinerant in the critical region [5]. Another pressure- induced 
superconductor among the Ce-Rh-Si ternaries is CeRhSis [1]. At ambient pressure, this 
tetragonal compound orders antiferromagnetically at Tn = 1.6 K and behaves as a Kondo 
lattice with the characteristic temperature Tk of about 100 K [6]. In contrast to the case 
of CeRh2Si2, in an applied pressure Tn does not change monotonically: it first increases 
up to 1.9 K near 0.8 GPa and then decreases down to 1.1 K near 2.6 GPa without 
further change at higher pressures [8]. Most importantly, superconductivity appears in 
this compound already under small pressures and coexists with the antiferromagnetic 
ordering in the entire pressure range studied, with the critical temperature gradually 
rising up to 1.1 K with increasing pressure [8]. Above 2.6 GPa, Tc merges with Tn 
[TJ [8] . Moreover, in the vicinity of this pressure, pronounced non- Fermi liquid features 
are observed in all the bulk properties. The occurrence of superconductivity in CeRhSis 
is highly peculiar, as the crystal structure of this compound lacks an inversion symmetry 
|8]. Apparently, CeRhSis is another member of a novel family of unconventional non- 
centrosymmetric superconductors with the most prominent representative being CePtsSi 
[9] . Likely, mixed spin-triplet pairing symmetry in CeRhSis is responsible for extremely 
large upper critical magnetic field -Bc2(0) = 30 T measured in a magnetic field applied 
along the crystallographic c-axis, its very unusual temperature dependence, as well as 
strong anisotropy of -Bc2(0) (cf. -Bc2(0) = 7 T for S ± c-axis) [TO] . 

Motivated by the intriguing physics in CeRh2Si2 and CeRhSis, we started 
systematic investigations of the physical behavior in other ternaries from the Ce-Rh- 
Si phase diagram. Recently, we reported on very complex magnetic properties of the 
orthorhombic compound CeRhsSi2, marked by two subsequent antiferromagnetic phase 
transitions at Tni = 4.7 K and Tn2 = 4.5 K, multiple metamagnetic transitions in 
the ordered state, and huge magnetocrystalline anisotropy being quite uncommon for 
Ce-based intermetallics [H]. In turn, for hexagonal Ce2RhSi3 we provided evidence of 
Kondo lattice behavior with the characteristic temperature Tk ~ 9 K, which coexists 
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below Tn = 4.5 K with long-range antiferromagnetic ordering [12]. Remarkably, by 
means of bulk and spectroscopic measurements, the cerium 4/ electrons were found to 
be well localized in both compounds. 

An entirely opposite character of the 4/ states may be expected for the silicides 
CeRhSi2 and Ce2Rh3Si5, which have been reported in the literature as valence 
fluctuating systems [131 EH ESI ES]. The former compound crystallizes with an 
orthorhombic crystal structure of the CeNiSi2-type. Its unit cell volume does not 
follow the lanthanide contraction established for the i?i?RhSi2 {RE = La, Pr, Nd) 
series. The magnetic susceptibility of CeRhSi2 shows a non- Curie- Weiss temperature 
variation with a broad maximum near 80 K. The electrical resistivity is proportional to 
the squared temperature below about 50 K, while above 100 K it decreases slightly with 
increasing temperature. These features of the bulk magnetic and electrical behaviors 
reflect instability of the Ce ions' valence. The Coqblin-Schrieffer approach applied to 
the magnetic data of CeRhSi2 yielded the characteristic temperature Tq = 309 K (Tq 
is related to the Kondo temperature Tk via the Wilson number, Tk = WTq). In turn, 
analysis of the resistivity data in terms of the Freimuth model gave estimates for the 
spin-fluctuation temperature Tgf = 147 K and the position of the 4/ band with respect 
to the Fermi energy Tf = 30 K. The valence fluctuating nature of CeRhSi2 has also been 
concluded in an independent study reported in Ref. [H]. The other compound of our 
interest, i.e. Ce2Rh3Si5, was studied in less details. It was reported to form with the 
orthorhombic U2Co3Si5-type structure as a member of the i?i?2Rh3Si5 {RE = La-Er) 
family. As the lattice parameters of its unit cell follow the lanthanide contraction along 
the series it was concluded in Ref. [15] that the Ce ions in this compound are trivalent. 
On the other hand, the same authors reported nearly temperature independent magnetic 
susceptibility of Ce2Rh3Si5, and in order to reconcile their contradicting flndings they 
speculated that the ground state in this compound is nonmagnetic because of the 
presence of 'virtual' spin fluctuations, which however do not destabilize the charge state 
of the Ce ions. The electrical transport properties of Ce2Rh3Si5 were reported in Ref. 
|16] . In the temperature range 2-30 K, the resistivity was shown to vary as T^, and 
interpreted as being governed by interband s — d scattering. Nonmagnetic character 
of Ce2Rh3Si5 was also inferred from the heat capacity data, which show featureless 
temperature behavior up to 30 K and strongly reduced entropy S'(30K) = 0.24R. 

In this work we reinvestigated the magnetic and electrical transport properties 
of CeRhSi2 and Ce2Rh3Si5. Moreover, we measured for the flrst time the speciflc 
heat of both materials in a wide temperature range 0.35-300 K, and performed 
x-ray absorption spectroscopy (XAS) and x-ray photoelectron spectroscopy (XPS) 
studies. The experimental data are discussed in the context of the calculated electronic 
structures. Altogether, the obtained results conflrm the valence fluctuating nature of 
both compounds. 
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2. Experimental and computational details 

Polycrystalline samples of CeRliSi2 and Ce2Rh3Si5 were prepared by arc melting 
stoichiometric amounts of the elemental components (Ce - 3N, Ames Laboratory, Rh 
- 3N, Chempur and Si - 6N, Chempur) in a copper-hearth furnace installed inside a 
glove-box filled with ultra-pure argon gas with continuously controlled partial pressures 
of O2 and H2O to be lower than 1 ppm. The buttons were flipped over and remelted 
several times to ensure good homogeneity. The weight losses after the final melting were 
negligible (less than 0.2%). 

The quality of the obtained alloys was checked by x-ray powder diffraction on an 
X'pert Pro PanAnalytical diffractometer with CuK^^ radiation and by energy dispersive 
x-ray (EDX) analysis using a Phillips 515 scanning electron microscope equipped with 
an ED AX PV 9800 spectrometer. Both techniques proved single-phase character of the 
Ce2Rh3Si5 sample, with the expected stoichiometry and crystal structure. In the case 
of CeRhSi2, however, some small admixture of CeRh2Si2 was evidenced in the x-ray 
pattern and the EDX spectrum. The structural refinements were done employing the 
program FULLPROF [T7] . 

Magnetic susceptibility measurements were performed in the temperature range 
1.72-800 K in magnetic fields of 0.5 T using a Quantum Design superconducting quantum 
interference device (SQUID) magnetometer. The heat capacity and the electrical 
resistivity were measured within the temperature interval 2-300 K using a Quantum 
Design PPMS platform. 

The x-ray photoelectron spectroscopy experiments were carried out at room tem- 
perature using a Physical Electronics PHI 5700/660 spectrometer with monochroma- 
tized AlKa radiation (1486.6 eV). The spectra were collected on parallelepiped-shaped 
specimens broken in-situ in high vacuum of the order of 10^° Torr. 

The measurements of the x-ray absorption spectroscopy at the Ce-Lni threshold 
were performed at several different temperatures at the EXAFS-1 beamline C of the 
Hamburger Synchrotronstrahlungslabor (HASYLAB/DESY) using a fixed-exit double- 
crystal Si (111) monochromator. In these studies Ce02 and CeFs were used as the 
internal standards. The energy scans were made step-by-step with an energy resolution 
of approx. 1 eV. The absorption spectra were calculated by ln{Ci/C2) from the x-ray 
intensities Ci, C2, detected by ionization chambers in front and behind the irradiated 
fiat sample. A peak to background ratio P/B = 0.6 at the Ce-Lni white line was 
realized with powdered samples of about 16 mg, which were ground with small amounts 
of B4C and fixed with paraffin wax on a 1 cm^ window of the fiat copper sample holder. 
Temperatures from the range 5-293 K were obtained by a He-gas flow cryostat showing 
thermal stabilities of T = -|-/ — 0.5 K during typical measuring dwell times of 20 minutes. 

The electronic structure calculations were performed within density functional 
theory [18] using the full-potential local-orbital minimum-basis band structure scheme 
(FPLO) [19] and the full potential linearized augmented plane wave (LAPW) method 
implemented in the latest version (WIEN2k) of the original WIEN code |20j. The 
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FPLO calculations were performed in the fuUy-relativistic mode and the LSDA (the 
local spin density approximation) exchange-correlation potential was assumed in the 
form proposed by Perdew and Wang [21]. In the Wien2k code calculations the 
scalar relativistic approach was implemented with the spin-orbit interactions taken into 
account using the second variational method [22]. Two different exchange-correlation 
potentials in the generalized gradient approximations (GGA) were tested, in the forms 
proposed by Perdew et al. [23j as well as Wu and Cohen [M]- Furthermore, to improve 
the description of the strongly correlated 4/ electrons, the on-site Coulomb energy U 
correction was introduced within the LSDA+U approach [25]. The value of Ug// was 
chosen equal to 6 eV [2^. The number of k-points was 8000 in the Brillouin zone 
(BZ), which corresponds to at least 1100 points in the irreducible wedge of the BZ for 
all systems and methods of calculations. For BZ integration a tetrahedron method was 
used [2l]. The self-consistent criterion was equal to at least 10~^ Ry for the total energy. 
The calculations were performed for lattice constants and atomic positions in the unit 
cells as given in the IIIA section. 

The theoretical photoemission spectra (XPS) were obtained from the calculated 
densities of electronic states (DOS) convoluted by Gaussian with a half-width (5) equal 
to 0.3 eV and scaled using the proper photoelectronic cross sections for partial states 



3. Results and discussion 

3.1. Crystal structures 

The compound CeRhSi2 was reported in the literature to crystallize with an 
orthorhombic structure of the CeNiSi2-type [131 [H]. In turn, Ce2Rh3Si5 was considered 
as having an orthorhombic structure isotypic to U2Co3Si5 [151 [16]. However, to the best 
of our knowledge, for none of these two phases were any details on the atomic positions 
in the unit cells reported. Therefore, the x-ray powder diffraction data obtained in the 
present study as part of characterization of the samples quality were used to refine the 
crystal structures of both compounds. 

For CeRhSi2, the refined lattice parameters are a = 4.2615(3) A, 6 = 16.7469(9) 
Aand c = 4.1751(3) A, in good agreement with those given in Ref. [H]. The crystal 
structure belonging the space group Cmcm was refined down to the residuals = 3.4% 
and i?„p = 5.3%. The obtained atomic coordinates are given in Table 1, together the 
values of equivalent isotropic thermal displacement parameters for all the atoms. The 
compound was thus confirmed to be isostructural with CeNiSi2. Detailed discussion of 
this crystal structure can be found in the literature, e.g. in the original paper [29]. For 
the purpose of this work it is enough to note that the unit cell contains one position of 
Ce atoms, which are coordinated by four Rh atoms at a distance of 3.226 A, four Sil 
atoms at a distance of 3.163 A, and two Si2 atoms at a distance of 3.167 A. 

As for Ce2Rh3Si5, the unit cell of the U2Co3Si5-type (space group Iham) has also 
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been corroborated. The crystal structure was refined down to the residuals -Rp = 2.3% 
and i?wp = 3.1%. The obtained lattice parameters are a = 9.8949(3) A,b= 11.7576(3) 
Aand c = 5.8114(1) A. These values differ only slightly from those reported in Refs. 
p3| [16] . The atomic coordinates and the equivalent isotropic thermal displacement 
parameters are collected in Table 2. In the unit cell of Ce2Rh3Si5 there is one position 
of Ce atom that is surrounded by one Rh2 atom at a distance of 3.138 A, two Rh2 atoms 
at a distance of 3.161 A, one Rh2 atom at a distance of 3.375 A, two Rhl atoms at a 
distance of 3.381 Aand one Rh2 atom at a distance of 3.399 A. The nearest neighbors 
silicon atoms are located at the distances: 3.022 A(one Si3 atom), 3.111 A(two Si3 
atoms), 3.129 A(one Sil atom), 3.148 A(one Si3 atom), 3.199 A(two Si2 atoms) and 
3.233 A(one Si2 atom). Further details on the U2Co3Si5-type crystal structure can be 
found in the literature (see e.g. Ref. [T6]). 



3.2. Magnetic properties 

The magnetic data of CeRhSi2 and Ce2Rh3Si5 are summarized in Fig. 1. The 
compounds exhibit small (of the order of 10~^ emu/mol per Ce atom) and weakly 
temperature dependent magnetic susceptibility that clearly manifests nonmagnetic 
character of both materials. With decreasing temperature from 800 K the susceptibility 
of CeRhSi2 slightly increases, goes through a broad maximum centered around 75 K 
and then rapidly rises below about 20 K. Similar behavior is observed for Ce2Rh3Si5, 
with a maximum located around 270 K and some tendency for saturation before the 
low-temperature upturn occurs. The overall shape of these x(T) curves is typical for 
Ce-based intermetallics with valence fluctuations. The position of the maximum in 
x(T), T(x™'"), gives an estimate for the characteristic temperature Tgf, related to spin 
fluctuations in such compounds. From the relation [30] 



one finds Tsf ~ 112 K for CeRhSi2 and Tgf ~ 405 K for Ce2Rh3Si5. Comparison 
of these values suggests that in these two compounds interactions of the cerium 4/ 
electrons with the conduction band is stronger in the latter one, in line with smaller 
and less temperature dependent magnetic susceptibility observed for Ce2Rh3Si5. 

Low-temperature upturns in x{T) curves are commonly observed in Ce-based 
intermediate valence materials and their origin is usually attributed to the presence 
of stable Ce^^ ions located at grain boundaries or/and some contamination by 
paramagnetic impurities. To account for this spurious effect it is assumed that the 
intrinsic susceptibility of CeRhSi2 and Ce2Rh3Si5 at low temperatures is given by the 
formula [31] 

X(0) = ^ (2) 

in which C = 0.807 emu/(mol K) stands for the Curie constant of free Ce^"*" ions 
{C = -^p^ where is the Avogadro number, ks is the Boltzmann constant, while /Xgff 
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= 2.54 is the effective magnetic moment of tlie cerium 4/^ state). The above-derived 
values of Tgf imply x(0) equal to 3.6 ■ 10~^ emu/mol for CeRhSi2 and 1.0 ■ 10~^ emu/mol 
per Ce atom for Ce2Rh3Si5. Then, the measured magnetic susceptibility can be modeled 
by the function 

X{T ^ 0) = x(0) + — (3) 

where the Curie term represents the contribution due to spurious Ce'^"*" ions in 
an amount of n atoms per mole of the given compound. Fitting this equation to the 
experimental data taken below 20 K (note the dashed lines in Fig. 1) yields n = 6.3-10"^ 
and 3.5-10"^ for CeRhSi2 and Ce2Rh3Si5, respectively. The magnetic susceptibility data 
corrected for the spurious Ce^"*" ions is shown in Fig. 1 by the full symbols. 

At higher temperatures the magnetic susceptibility of both compounds can be 
analyzed in terms of the interconfiguration fluctuations (ICF) model developed for 
intermediate valence systems by Sales and Wohlleben [32]. Within this approach the 
magnetic susceptibility of a Ce-based compound with nonmagnetic Af^ ground state is 
expressed as 

J + Jo 

where Tq is a characteristic temperature associated with valence fluctuations 
between the 4/° and 4/^ configurations of Ce ions, whereas z^(T) stands for the 
temperature dependent mean occupation of the ground state that is given by the formula 

1 + 6exp(^) 

in which i?ex denotes the energy difference between the ground and excited states. 
The constant term Xo iii Eq. 4 accounts for paramagnetic contribution due to conduction 
electrons and diamagnetic contribution due to core electrons. Fitting the ICF model to 
the experimental data of CeRhSi2 and Ce2Rh3Si5 resulted in reasonable descriptions of 
the magnetic behavior above about 50 K and 150 K, respectively (note the solid lines 
in Fig. 1). The parameters obtained in the least-squares fits are Eq^ = 220 K, Tq = 
44 K and Xo = 2.3 • 10~^ emu/mol for the former compound, and i?ex = 845 K, Tq 
= 129 K and xo = 4.9 • 10~^ emu/mol for the latter one. Assuming that the highest 
possible valence for the cerium 4/° state is +3.30 (for a discussion see Ref. [33]), one 
may estimate from Eq. 5 the change of the effective valence of Ce ions with varying 
temperature. In the range in which the ICF model provides proper approximation of 
the susceptibility of the compounds studied one finds for CeRhSi2 the valence change 
from +3.19 at 50 K to +3.07 at 300 K and to +3.05 at 800 K, while for Ce2Rh3Si5 the 
change is from +3.23 at 150 K to +3.16 at 300 K and to +3.09 at 800 K. 
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3.3. Heat capacity 

Fig. 2 displays the temperature dependences of the specific heat of CeRhSi2 and 
Ce2Rh3Si5. The sohd hnes through the experimental points are the least-squares fits of 
the formula: 

C(T) = ai(T) + Cph(T), (6) 

where the first term is the electron specific heat Cei(T) = 7T, while the second one 
represents the phonon contribution in the form [MJ |35] : 



CAT) 



1 - aT 



V— 

ph,E 



Here, a stands for the anharmonic coefficient, R is the gas constant, Cph,D describes 
the Debye contribution of three acoustic modes (characterized by the Debye temperature 
9d), and Cph,E represents the Einstein specific heat calculated for i groups with Ui 
optical branches (characterized by the Einstein temperatures Oej)- In order to avoid 
overparametrization of the experimental curves C{T) (Fig. 2), we reduced the number 
of the groups of optical branches to z = 2, which seems to be minimum to describe 
satisfactorily the specific heat data. In that case we found that the best results one 
obtains assuming multiplicity rii = 2 and n2 = 7 in CeRhSi2, and ni = 12 and n2 = 15 
in Ce2Rh3Si5. The so-obtained explicit form of Eq. (E]) applied to the experimental 
data above about 15 K yields for CeRhSi2 the fitting parameters: 6d = 163 K, 9ei 
= 139 K, = 324 K, 7 = 22.4 mJ/(mol K^) and a ^ 1.0 ■ 10'^ 1/K. In the case 
of Ce2Rh3Si5 a reasonable fit to the experimental data was obtained with 0d = 147 
K, Gei = 198 K, Ge^ = 427 K, 7 = 45.3 mJ/(mol K^) and a = 2.7 • 10"^ 1/K. It is 
worthwhile emphasizing that the presented model does not reflect complexity of the real 
lattice vibrations in the two compounds, and at most gives only basic characteristics of 
the phonon spectra. 

As can be inferred from the insets to Fig. 2, below about 12 K the experimental 
specific heat data can be described by a simple formula: 

1944 

el 



C(T)/T = 7 + r-^T^ (8) 



'D,LT 

where the second term is the conventional T^-Debye law with the low-T Debye 
temperature 0d,lt and r stands for a number of atoms in a formula unit [M| [35] . The 
insets to Fig. 2 present the least-squares fits of the latter formula to the measured 
specific heat data with the parameters 7 = 86 mJ/(mol K^) and 6d,lt = 315 K for 
CeRhSi2 and 7 = 47 mJ/(mol K^) and 6d,lt = 375 K for Ce2Rh3Si5. The obtained 
values of Od,lt are obviously different from those obtained using previous full-range 
fittings since in the latter approximation only the acoustic branches are considered, and 
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thus 0D,LT reflects an average characteristic temperature of the phonon spectrum. The 
obtained 7 coefficients are in both compounds of the order typical for systems with 
valence fluctuations [36]. 

At low temperatures approaching the experimental limit a little upturn is observed 
in C{T)/T of CeRhSi2. This anomaly has an unknown origin and can temporarily be 
ascribed to the impurity phase CeRh2Si2, detected on the x-ray diffraction pattern of the 
sample measured. However, an intrinsic nature of the upturn can neither be excluded. 

3.4- Electrical resistivity 

The temperature variations of the electrical resistivity of CeRhSi2 and Ce2Rh3Si5 are 
shown in Fig. 3. For both compounds the absolute magnitude of the resistivity is 
very high because of many cracks present in the specimens measured (the samples 
were very brittle). Therefore, quantitative discussion of the electrical behavior of these 
silicides is not possible. Nevertheless, it is clear that both materials exhibit metallic- 
like conductivity with some features characteristic of valence fluctuations systems [36] . 
Below 50 K for CeRhSi2 and 90 K for Ce2Rh3Si5 K the resistivity is proportional to 
the squared temperature, as predicted for such systems within the paramagnon model 
[30] (some deviations from p ~ are observed below 10 K and 20 K, respectively). At 
higher temperatures the p(T) curve of CeRhSi2 forms a broad shallow maximum above 
100 K, while that of Ce2Rh3Si5 shows a tendency to saturate near room temperature. 
The observed behavior is thus fully consistent with the values of the spin fluctuation 
temperature T^f ^ 112 K and Tgf ~ 405 K estimated from the magnetic susceptibility 
data for the former and the latter compound, respectively. 

3.5. XAS spectra 

The x-ray absorption spectroscopy data obtained for CeRhSi2 and Ce2Rh3Si5 are 
gathered in Fig. 4. For each compound a linear baseline determined at 100 eV below 
the Ce-Liii absorption edge was subtracted, and the spectrum was normalized to the 
intensity values at E = 5755 eV. The absorption spectra of both phases show one 
dominant white line at 5726 eV of the Ce-Lni line and several maxima with energies 
higher than 5750 eV in the extended x-ray absorption fine structure (EXAFS) region. 
The EXAFS signal is influenced by specific features of the crystal structure of the 
investigated compound. Comparison of the measured spectrum with that of the Ce^^ 
reference system CeP04 allows us to attribute the observed white line to the same 
dominant 4/^ electronic configuration of the cerium atoms. Both phases, i.e. CeRhSi2 
and Ce2Rh3Si5, thus show a mean valence of cerium ions relatively close to +3. However, 
the presence of an additional contribution at nearly AE ^ 9 eV above the white 
line maximum clearly shows mixed valent behavior of the compounds investigated. 
Decreasing the temperature from room temperature to T = 5 K leads to a small decrease 
in the intensity of the white line and simultaneous increase of the other peak. The 
observed shift of the spectral weight of the Ce-Lm white line to higher energies implies a 
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reduced influence of tlie 4/ electron on the 2p — ?■ nd transition. This behavior manifests 
a gradual gain in the contribution due to the 4/° configuration, and thus reflects an 
increase of the mean valence of the Ce atoms with decreasing temperature. 

In quantitative evaluation of the spectra measured for CeRhSi2 and Ce2Rh3Si5, 
two Gaussian functions were considered, which represent the 4/^ and 4/° contributions, 
respectively (see an example presented in panel (c) of Fig. 4). Additionally, an arctan- 
function was added in order to account for the contribution due to transitions of the 2p 
electrons to the conduction band. Analyzing the relative increase of the 4/° peak with 
decreasing temperature one obtains for CeRhSi2 the mean valence of Ce ions changing 
from +3.15 at room temperature to +3.17 at 5 K (see Fig. 4d). Similar calculation 
performed for Ce2Rh3Si5 yields the valence increase from +3.18 at 280 K to +3.22 at 5 
K (cf. Fig. 4b). 

3.6. XPS spectra 

The XPS spectra of the 2id core levels usually provide some detailed information 
about the 4/ shell configurations and the /-conduction-electron hybridization. The 
calculations by Gunnarsson and Schonhammer [38] indicate that a large hybridization 
potential, V/s, is necessary to explain the observed spectra for several Ce-intermetalhc 
compounds, which often show different final states depending on the occupation of the 
/ shell: /°, and p (see Ref. [39]). Due to the spin-orbit (SO) interaction there are 
two sets of Ce ?>d photoemission lines in the spectrum attributed to the 3^3/2 and 3^5/2 
components of the final states, with an intensity ratio /(3(i5/2)//(3(i3/2) = 3/2. The 
main photoemission lines originating from Ce^"*" are labeled as 3(i4/^. 

The SdAf"^ final state components appears when the core hole becomes screened by 
an additional 4/ electron, which is possible due to the hybridization of the Ce 4/ shell 
with the conduction electrons. Consequently, the 3c?^/^ components in the Ce 3d XPS 
spectra are attributed within the Gunnarsson-Schonhammer model to the /-conduction 
electron hybridization energy Afg. Afs = ttV^ pmax describes the hybridization part of 
the Anderson impurity Hamiltonian [IQ], where Pmax is the maximum in the DOS and 
V is the hybridization matrix element. It is possible to estimate Afs from the ratio 
r = /(/^)/(/(/^) + J(/^)), calculated as a function of Afs in Ref. [39], when the peaks 
of the Ce 3d XPS spectra that overlap are separated. 

Fig. 5 shows the Ce 3d XPS spectra of CeRhSi2 (a) and Ce2Rh3Si5 (6). The 
separation of the overlapping peaks in the spectra was done on the basis of Doniach- 
Sunjic theory [H]. A background, calculated using the Tougaard algorithm [12], was 
subtracted from the XPS data. The estimated value of the SO splitting equals 18.6 
eV. In Fig. 5 each SO set of the Ce 3d photoemission lines consists of contributions 
marked as /°, and The intensity ratio r ^ 0.09 for CeRhSi2 and ~ 0.17 for 
Ce2Rh3Si5. This intensity ratio gives for CeRhSi2 and Ce2Rh3Si5 a crude estimate of a 
hybridization width Afs ~ 40 meV and ~ 83 meV, respectively. The estimated energy 
Afs for Ce2Rh3Si5 is twice as large as the value for CeRhSi2, which suggests a lower 
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occupation number of the 4/ shell in Ce2Rh3Si5. 

In order to determine the ground state / occupation from the ?>d XPS spectra we use 
Fig. 4 and 6 of Ref. [39], where the dependence of the ratio /[{I{f) + I{p) + I{p)] 
on the / occupation is shown for different Afg. Relative /° intensities of magnitude 
~ 0.05 for CeRhSi2 and ~ 0.07 for Ce2Rh3Si5 corresponds to n/ value of the order 
0.95 and 0.93, respectively. The Ce ?>d XPS spectra, however, allow an estimate of the 
occupation number Uf (and of the energy Afs) within an accuracy of the order of 20%. 
The errors due to the uncertainties in the intensity ratios we discussed previously [13] . 

We note, however, that the Uf value obtained from the Ce ?)d XPS spectra is 
smaller than that estimated from the X-ray absorption spectra either for CeRhSi2 or 
Ce2Rh3Si5. This discrepancy show immediately that the final-state / occupation is a 
different function of the initial-state occupancy for each type of spectroscopy (see Ref. 
^44j ) . We note, however, that the Uf value for CeRhSi2 is larger than uj derived for 
Ce2Rh3Si5. 

Other evidence for an fluctuating valence in Ce2Rh3Si5 comes from the Ce Ad XPS 
data. As shown in Fig. 6, the latter spectrum exhibits some features above 120 eV 
(marked by the arrow), which can be assigned to the /° final state [l5l 116] . 

3. 7. Electronic hand structure 

The band structure calculations performed using two distinct approaches, i.e. with 
FPLO and FP-LAPW, yielded very similar results and therefore only those obtained by 
using the Wien2k code are discussed below. The latter method gives also an opportunity 
to go beyond LDA using L(S)DA-|-U approach. The calculations were performed with 
and without spin polarization. The spin-polarized calculations were started with finite 
initial magnetic moments in order to promote magnetic solutions. However, the self- 
consistent results for magnetic solutions were nearly same as those for nonmagnetic 
ones. In all the tests performed, the magnetic moments on cerium atoms were always 
below 10^^ yUs/Ce atom for both CeRhSi2 and Ce2Rh3Si5. 

The calculated densities of states (DOS) for the compounds CeRhSi2 and Ce2Rh3Si5 
are presented in Figs. 7 and 8, respectively. In both cases the DOS plots can be divided 
into three parts below the Fermi level {Ep = 0). The first part located below -16 eV is 
composed of two peaks formed by the Ce(5pi/2) and Ce(5p3/2) electrons. It should be 
noticed that the LDA+U peaks are shifted towards higher binding energies. As will be 
discussed below, this shift improves consistency between the calculated and measured 
X-ray photoelectron spectra, hence justifying the LDA+U approach. The second part, 
located between about -6 and -11 eV, is formed mainly by s electrons of the Rh and Si 
atoms. The third part, which is the main one of the valence bands, is located between the 
Fermi level and about -6 eV of binding energy. The main contribution to these subbands 
is provided by d electrons located on the Rh atoms and Si(3p) electrons. A small 
contribution comes also from Ce(5d) electrons, but less than one electron per Ce atom. 
The Ce(4f ) electrons form a narrow band with the center of gravity above the Fermi level. 
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Most of the band is empty, and within the LDA scheme the number of occupied Ce(4f) 
states below is equal to 0.92 for both systems. The LDA+U calculations change 
the electronic structure in a way that the Ce(4f) band is slightly shifted above Ei? and 
consequently the number of Ce(4f) electrons is reduced to 0.57 and 0.61 for CeRhSi2 
and Ce2Rh3Si5, respectively. For the former compound, the occupation of electrons 
inside the atomic sphere of Ce atom is given by the LDA approach as 6s^'^^5d*^'^^4f°'^^ 
and by the LDA+U calculations as 5si-995cl0-8i4f0-57_ xhese results should be compared 
with the configurations 6s^5d^4f^ and 6s^5d^4f^ of Ce"*"^ and Ce"*"^ ions, respectively. 
Apparently, the calculated occupancies are intermediate between those characteristic 
for stable valence ions. Similarly, the results were obtained for Ce2Rh3Si5, which are 
ggi.985^o.674f0.92 Qgi.995jo.824fO.6i and LDA+U calculations, respectively, 

also indicate distinct deviation from the stable valence configuration. This effect may 
be caused by additional charge transfer from/to rhodium and silicon atoms as well as 
by the charge accumulated in interstitial regions between atomic spheres. 

The total DOS at Ei? and the main contributions of particular atoms are 
summarized in Table 3. The main contributions, about 50%, to the total DOS at 
Ei? are provided by Ce(4f) electrons. Based on the DOS at E^? one can calculate 
the Sommerfeld coefficient 70, to be equal to about 11.6 and 11.2 mJ/(mol K^) for 
CeRhSi2 and Ce2Rh3Si5, respectively, obtained within the LDA approach and 3.7 and 
6.5 mJ/(mol K^), respectively, obtained with the LDA+U approach. Comparison 
of these values to the experimental ones (7 = 86 and 47 mJ/(mol K^, for the two 
compounds respectively) yields the mass enhancement factors A=(7exp/7o)-l, which are 
also collected in Table 3. 

Fig. 9 presents the X-ray photoelectron spectra of CeRhSi2 and Ce2Rh3Si5, 
calculated from the DOS data using the cross sections reported in Ref. p8j. For 
comparison, the experimental XPS spectra are also shown. Apparently, the calculated 
results are quite consistent with the experimental ones. It should be noted that the 
LDA+U approach better describes the measured spectra for both systems, as judged 
from the distinct differences in the region of the Ce(5p) peaks. 

4. Summary 

The macroscopic and spectroscopic results obtained in this work for CeRhSi2 
and Ce2Rh3Si5 unambiguously corroborate intermediate valent character of both 
compounds. Mutually comparing the properties of these two phases it is worth noting 
that the effective valence of Ce ion in CeRhSi2 is closer to +3 than in the other silicide 
in the entire temperature range studied. In both compounds the electronic ground state 
of cerium is 4/°, yet the energy distance to the excited magnetic 4/^ state, £'ex, is 
smaller in CeRhSi2 and hence the magnetic susceptibility of this compound is larger 
and more temperature dependent than the susceptibility of the other one. Furthermore, 
stronger hybridization of the 4/ electronic states with the conduction states in Ce2Rh3Si5 
is clearly reflected in larger values of the characteristic spin and charge fluctuation 
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temperatures Tq and Tgf. Similar conclusions as regards the valence states in the two 
compounds can be derived from the XAS and XPS results. The mean occupation of 
the 4/ state is smaller in Ce2Rh3Si5 and less temperature dependent than in the other 
silicide, in line with the bulk magnetic data. The experimental XPS spectra can be 
fairly well reproduced by ab-initio electronic band structure calculations, which yielded 
nonmagnetic ground state in both compounds with significant mass enhancements due 
to strong electronic correlations. 
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Table 1. Atomic positions and equivalent isotropic thermal displacement parameters 
for CeRhSia. 

y z U(eq) 



Atom 


Site 


X 


Ce 


4c 





Rh 


4c 


0.5 


Sil 


4c 





Si2 


4c 


0.5 



0.8945(1) 0.75 0.0108(6) 

0.8211(1) 0.25 0.037(1) 

0.5427(7) 0.75 0.032(4) 

0.7477(5) 0.75 0.010(3) 



Table 2. Atomic positions and equivalent isotropic thermal displacement parameters 
for CeaRhaSis. 



Atom 


Site 


X 


y 


z 


U(eq) 


Ce 


8j 


0.2668(1) 


0.1305(1) 





0.0072(3) 


Rhl 


4a 








0.75 


0.0097(6) 


Rh2 


8j 


0.3923(1) 


0.1386(2) 


0.5 


0.0108(4) 


Sil 


4b 


0.5 





0.75 


0.018(2) 


Si2 


8g 





0.2218(4) 


0.75 


0.017(2) 


Si3 


8j 


0.1574(6) 


0.1095(6) 


0.5 


0.021(2) 



Table 3. DOS at E^: total (per eV and f.u.) and site projected (per eV and atom) 
from orbitals giving the main contributions to the total value; calculated Sommerfeld 
coefficients 70 in mJ/(mol K^) and mass enhancement factors A. 



Compound 


Type od 


Total 


7o 


A 


Ce(4f) 


Rhl(4d) 


Rh2(4d) 




calculations 














CeRhSia 


LDA 


4.92 


11.60 


6.4 


2.87 


0.32 






LDA+U 


1.58 


3.72 


22.1 


0.67 


0.16 




Ce2Rh3Si5 


LDA 


4.73 


11.15 


3.2 


1.19 


0.55 


0.19 




LDA+U 


2.74 


6.46 


6.3 


0.61 


0.27 


0.15 
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Figure 1. Temperature dependences of the molar magnetic susceptibility of CeRhSi2 
and Ce2Rh3Si5. The full symbols represent the x(T) curves corrected for contributions 
coming from spurious Ce'^+ ions (see text for details). The dashed and solid lines are 
the fits discussed in the text. 



120 




Temperature (K) 

Figure 2. Temperature dependences of the specific heat of (a) CeRhSi2 and (b) 
Ce2Rh3Si5. The insets present the low-temperature data in the form C/T vs. . The 
solid lines are the fits discussed in the text. 
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Figure 3. Temperature dependences of the electrical resistivity of CeRliSi2 and 
Ce2Rh3Si5. The solid lines are the fits discussed in the text. 
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Figure 4. X-ray absorption spectra at the Ce-Lm threshold measured at various 
temperatures for (a) CeRhSi2 and (b) Ce2Rh3Si5. The standard spectrum for Ce'^"'" 
ions, i.e. that of CeP04 is represented by the dashed curves. As an example, 
deconvolution of the spectrum measured for CeRhSi2 at 5 K is shown in panel (c). 
The contributions due to 4/^ and 4/° configurations are depicted by the solid lines, 
while the dashed line represents the arctan-function accounting for transitions of the 
2p electrons to the conduction band. Panel (d) displays the calculated changes of the 
Ce ions valence in both compounds with varying the temperature. 
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Figure 6. Ce 4d XPS spectrum obtained for Ce2Rh3Si5. 
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Figure 7. The total (per f.u.) and site projected (per atom) densities of electronic 
states (DOS) for CeRhSi2 calculated within the LDA+U scheme (for comparison the 
total DOS and the Ce-atom DOS calculated using LDA approach are presented by the 
broken line). 
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Figure 8. The total (per f.u.) and site projected (per atom) densities of electronic 
states (DOS) for Ce2R3hSi5 calculated within the LDA-I-U scheme (for comparison the 
total DOS and the Ce-atom DOS calculated using LDA approach are presented by the 
broken Hne). 
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Figure 9. XPS valence band spectrum of (a) CeRhSi2 and (b) Ce2Rh3Si5, calculated 
within the LDA and LDA+U approaches and compared to the experimental one 
(represented by the symbols). Note the difference between the calculated spectra in 
the region of 5p contribution. 



